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ABSTRACT 

A mathematical model has teen developed for the 
absorption of carbon dioxi de from ammonia synthesis gas in 
amine promoted hot potassium carbonate solution. The model 
uses vapour— liquid— equilibrium data for the unpromoted 
carbonate~CC> 2 system which seems to be adequate at low 
fractional saturation of the liquid with carbon dioxide. In 
the temperature range of industrial operation, instantaneous 
reaction rate provides the appropriate enhancement factor for 
the chemical absorption. The calculated results have been 
found to be in close agreement with the observed values 
for an industrial unit. The model has been used to calculate 
various flow, temperature and concentrat ion profiles in the 
absorber. Variation of exit C0 2 concentration with seme of 
the operating parameters has also been examined. 



CHAPTER 1 


INTRODUCTION 

For a commercially useful ammonia synthesis process, 
economic removal of carbon dioxide from the synthesis gas 
has considerable importance. Carbon dioxide acts as a 
diluent of the synthesis gas and is also thought to be a 
poison for the synthesis catalyst when present in large 
concentration. Industrially, removal of carbon dioxide is 
largely brought about in two steps. The first step uses an 
absorption system where a chemical solvent is employed to 
reduce C0 2 to about 0.1% by volume. In the second step, a 
catalytic gas-phase methanation reduces total carbon oxides 
(COj + CO) to less than 10 ppm, as required to protect the 
ammonia synthesis catalyst. The first step of the sequence is 
responsible for bulk removal of CC> 2 and hence forms an important 
part of ammonia synthesis plant. The chemical solvents usually 
employed in industry are either the amines (e.g., mono—, di>», 
and triethanolamines) or hot potassium carbonate solution, 
usually promoted by organic amines. Though amine systems have 
a greater capacity for carbon dioxide and lower consumption 
of regeneration steam than the original alkali carbonate 
process operating at atmospheric pressure, they do not offer 
any significant advantage at elevated pressures, like the 
"hot potash" system (Benson et al., 1954), Moreover, hot 
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potash systems have higher energy efficiency owing to the 
fact that the heat exchange between the hot and cold streams 
can be eliminated. These have made the hot potash system 
evidently popular over the amine systems. Among the processes 
employing hot potassium carbonate as the basic chemical solvent , 
namely, the Benfield process, the Giammarco— Vetrocoke process, 
the Catacarb process and so on (Slack and James, 1974; Kohl 
and Risenfield, 1979; Strelzoff, 1981) the former has gained 
the broadest industrial acceptance for acid gas removal. The 
original Benfield process, developed at the U.S. Bureau of Mines 
in the early 1950s (Benson et al., 1954, 1956; Field et al,, 

1962) by Benson, Field and co-workers with a view to economically: 
remove CO^ and at high pressure frcm coal derived synthesis 

gas, has undergone several major modifications (Parrish and 
Ne±lson,l974; Bartoo and Ruzicka, 1983; Bartoo, 1984), all 
geared towards decreasing the investment costs and the utilities 
requirements while, at the same time, enhancing operational 
reliability and flexibility. At present, there are over 520 
installations throughout the world using the Benfield process 
with various modifications (Bartoo, 1984) and this indeed 
speaks for itself, the need to model the Benfield absorption 
system, more than any other system, as an integral part of an 
Ammonia plant Flowsheeting programme. A separate, yet concise, 
discussion on the Eenfield system, relevant to present work, 
is therefore given in Chapter 2. 
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The most important aspect of the Benfield system., that 
discerns it from other absorption systems is that the absorption 
of carbon dioxide is accompanied and enhanced by simultaneous 
chemical reaction with potassium carbonate. This makes the 
problem of modelling the system quite involved and prone to 
cumbersome calculations. Though the theory of chemical 
absorption was established in the late 60s (Danckwerts .. & Sharma, 
1966; Astarita, 1967; Ramm, 1968) adequate vapour-liquid 
equilibrium (VIE) data and the kinetics of reaction with con- 
centrated hot potash solution, particularly, for the high 
temperature range usually employed by the industry, were still 
not to be found in the open literature. Most of the data were 
proprietary and well beyond the reach of the scientists in the 
open world. Moreover, the effect of various organic promoters 
on the kinetics and VIE curve were not known in its entirety. 

All these had contributed to making the modelling of Benfield 
absorber a difficult, though not impossible, proposition. In 
early 70s, Uchida et al. (1972) published a rather simplistic 
"compartment-in— series" model for design and simulation of an 
adiabatic packed absorber, which was tested for the case of 
CO 2 absorption in hot potash solution. The thermodynamics 
and kinetics developed for the problem were anything but 
rigorous. Apart from this a computer program for design of 
hot carbonate absorption system was published by (Wellman and 
Kate 11 (1968). No mathematical model, however, was presented 
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aiongwith the package. Availability of recent accurate 
thermodynamic (Kent and Eisenberg, 1975; Astarita et al., 1983) 
and kinetic data (Astarita et al., 1983) as well as knowledge 
of the effects of rate promotion (Astarita et al., 1981; 1983) 
and other information related to this system, has paved the 
way for attempting simulation of the Benfield absorbers. 

The mathematical modelling of the Benfield absorbers 
with the relevant thermodynamic and kinetic models has been 
described in Chapter 3. The packed absorber has been modelled 
on the basis of two-film theory of mass transfer which for all 
practical purposes can be deemed at least as good as the more 
sophisticated penetration or surface renewal models (Astarita, 
1967). Other relevant parts of the model have been developed 
primarily on the basis of the theoretical analysis of Astarita 
et al. (1983). Recently published thermodynamics data and 
transpqrt properties have been incorporated in the present model 
in the form of suitable correlations obtained either directly 
fran literature or by means of regression analysis of available 
curves. The two— point boundary value problem., generated from, 
the steady state heat and material balance around a differential 
section of the packed absorber has been solved by the shooting 
method (Pox, 1957). The model developed in the present work 
has been validated with plant data. The fourth chapter 
presents a discussion of the results and gives the profiles 
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of flow rates, temperature and composition along the 
tower height for the gas and liquid streams. The implications 
of these and the scope of future modifications have been 
included in Chapter 5. A User's Manual for the computer 
program has been provided as an Appendix b. Program listing 
is given in Appendix C. 



CHAPTER 2 


THE BENFIELD PROCESS 

. s' ( 

2.1 Process Description 

Historically, the development of the famous "Benfield" 
process, named after its discoverers Benson and Field of U.S. 
Bureau of Mines, was stimulated in the early 1950s, by the 
urge of finding an improved method of removing carbon dioxide 
from Fischer—' Tropsch synthesis gas which would reduce the 
high capital and operating costs (Field et al., 1962). Since, 
coal, oxygen and steam are reacted under high pressure, at 
about 20 to 35 atm. in that process, it is more economical to 
use the absorption system at elevated pressure than at the 
atmospheric to eliminate a subsequent compression system for 
gases. This' holds equally good for ammonia synthesis gas which 
reacts under high pressure (15 to 25 atm) in the synthesis 
reactor. In order to exploit the advantage of the elevated 
pressure and temperature at which the impure gases are available, 
the use of hot concentrated alkaline solutions under pressure 
was investigated by the U.S. Bureau of Mines. In this system, 
the driving force which is the difference between the partial 
pressure of carbon dioxide in the bulk of the ,gas and its 
equilibrium partial pressure at the gas-liquid interface, was 
believed to be sufficient for absorption without the conven- 
tional lowering of temperature. Contrary to the usual practice. 
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absorption would then be possible at about the same temperature 
as regeneration. Since the carbonate-bicarbonate solution 
remains hot at all times, more concentrated solutions can be 
used than with cold absorption without precipitating the 
bicarbonate (Benson et al., 1954; 1956). 

Since a high concentration at alkalinity was desirable, 
potassium carbonate was considered to be more suitable than 
the sodium salt because the solubility of potassium bicarbonate 
was much greater than that of sodium bicarbonate. The bicarbo- 
nate is formed during absorption according to the following 
overall reaction 

K 2 C0 3 + C0 2 + H 2 0 TZ± 2KHC0 3 (2.1) 

Moreover, at normal temperature, the equilibrium partial 
pressure of carbon dioxide over solution of potassium carbonate 
was lower than that for comparable concentration _of sodium salt 
solution. Hence potassium carbonate was again preferred to 
sodium carbonate (Harte, et al., 1933). Figure 2.1 shows a 
simplified flowsheet for the potassium carbonate process. The 
absorber operates at an elevated pressure of about 25 atm and 
the regenerator at about 1 atm. The lean carbonate solution 
flows from the bottom of the regenerator and is pumped, without 
cooling, to the top of the absorber. The rich solution 
leaving the bottom of the absorber flows to the top of the 




PURIFIED 



Fig. 2.1 Simplified flow-sheet of single -stream absorber 
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regenerator afters pressure let, down. Seme dissolved ear bon 
dioxide and steam are flashed from the depressurized solution, 
but the rest of the carbon dioxide is stripped by steam 
produced by boiling the solution in the reboiler. The typical 
temperature of the rich solution leaving the absorber ranges 
from about 105 to 120°C, whereas the lean solution is at about 
120-125°C. In the scheme shown in Figure 2 «L, only single stream 
flows both in the absorber and the regenerator. This is 
generally used for carbon dioxide removal to a typical purity 
of 1% C0 2 . 

For a higher degree of purification, the flowscheme 
shown in Figure 2.2(a), a "Split-flow absorber" design is 
used. The bulk of the hot, regenerated lean solution, 60% - 
70% of the total, is introduced into the absorber at an 
intermediate point, while a small trim solution is cooled by 
10°'— 40°C and pumped to the top of the teswer. The lower tempera- 
ture trim stream has a reduced equilibrium back pressure of 
carbon dioxide to permit attaining high purity of 0.1% CO^ in 
the purified gas. Since the bulk of the absorber still 
operates at high temperature, the rich solution stream is hot 
and regeneration efficiency is not impaired (Bartoo, 1984). 

When still higher purity is required, a "two-stage" 
design, i.e., a split— stream flow in the absorber and split- 
stream flow in the regenerator is required (Figure 2.2(b)). 




Fig. 2.2 (a) Simplified flowsheet of cooled split-flow 
absorber 

2.2 (b) Simplified flowsheet of split- flow absorber 
with two stage regenerator 
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Two regenerated streams leave the regenerator; the major 
portion (60—85% of the total) is removed from an upper level 
in a partially regenerated state for return to an intermediate 
level of the absorber* The minor stream, 15—40% of the total 
is regenerated further in the bottom section. ..of the regenerator 
by the total boil, up steam, removed from the sump, and cooled 
before being pumped to the top of the absorber. In this scheme, 
the CO 2 level in the purified gas is brought down to 0.05%. 

The CO 2 cleanup attainable with the various flowsheet confi- 
gurations described above is summarized in Table 2.1. 

Table 2.1: CO 2 Cleanup Attainable with Promoted 


Hot Carbonate 

Process(*) 

Configuration 

Final CC >2 content 

Single stream absorber 

1% or greater 

Cooled split flew absorber 

0.1 - 1% 

Cooled split flow absorber 

and two— stage regenerator 

50CW1.000 ppm 


(*) A’starita, Savage & Bisio; Gas Treating with 
Chemical Solvents, 1983 


Recently other improvements over these conventional 
schemes have been proposed (Parrish and Neilson,l974; Bartoo 
and Ruzicka, 1983; Bartoo, 1984) fop example, the Hipure & 
Loheat Processes. The Hipure process employs two independent 
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solution circuits, arranged to have multiple use of the 
regeneration steam from the secondary circuit. The Loheat 
process recovers heat internally by reducing pressure over 
the regenerated lean solution and compressing the steam 
flashed off. 

For further details about various gas treating processes, 
the reader is advised to refer to "Ammonia" by Slack and James 
(1974), "Gas Purification" by Kohl and Risenfield (1979), 
"Technology and Manufacture of Ammonia" by Strelzoff (1981), and 
"Gas Treating with Chemical Solvent" by Astarita et al, (1983). 

2.2 Chemistrv of the Process: 

When carbon dioxide is absorbed in an unpromoted solution 
of potassium carbonate and bicarbonate, it reacts according to 
the following overall reaction 

C0 2 + K 2 C0 3 + H 2 0 — ■> 2KHCC> 3 (2.2) 

Since potassium carbonate and bicarbonate are both strong 
electrolytes, it may be assumed that the metal is present only 
in the form of K + ions, so that the reaction (2.2) may be more 
realistically represented in the ionic terms as 

C0 2 + C0 3 + W 2 C> 2HCC> 3 (2,3) 

In absorption, reaction 2.3 proceeds from left to right. The 
above reaction, being trimolecular ,is evidently made up of a 
sequence of elementary steps. Carbon dioxide may undergo 
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two direct reactions in an alkaline solution/ namely 

C0 2 + H 2 0— =?HCO^ + H + (2.4) 

C0 2 + OH 1 ^=£5' HCO 2 (2.5) 

Both reactions (2.4) and (2.5) are followed by an instantaneous 
reaction leading to the overall result developed by reaction (2.3). 

co^ + H + ^rr ± hc °3 (2.6) 

CO^ + H„0 V HCO^ + OH 1 (2.7) 

Both these reactions are accompanied with the following 
instantaneous dissociation reaction 

H 2 0 -a H + + OH 1 (2.8) 

The first sequence, namely, reactions (2.4), (2.6) and (2.8) 
is known as the "acidic" mechanism by Wall (Savage et al., 1980), 
The acidic mechanism's contribution to the overall reaction 
rate is negligible unless the pH of the liquid solution is 
very low. Almost all the cases of industrial absorption are 
held at high pH (pH >8) and hence the acidic mechanism is 
neglected altogether in the present study. Since reactions (2.7) 
and (2.8) are instantaneous reactions, reaction( 2.5 ) is the 
rate controlling step for absorption of C0 2 in unpromoted hot 


potash solution. 
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.2,3 Chemistry of Rate Promotion: 

Though significant rate enhancement is observed for 
the K^CO^ system of temperatures exceeding 318 K, the rate 
of chemical reaction is still not large enough, even at 
temperatures as high as 378 K, for the reaction to be regarded 
as instantaneous (Savage et al., 1980), Therefore, even at the 
highest temperatures of industrial practice the use of 
additives for rate promotion is of significant interest. 

Promotion of mass transfer rates, through additives to 
the K^CO^ solution, has been reported in the literature (Shrier 
and Danckwerts, 1969; Roberts and Danckwerts,. 1962; Danckwerts 
and Sharma, 1966; McLachlan and Danckwerts, 1972; McNeil and 
Danckwerts, 1964, 1967; Jeffreys and Bull, 1964; Astarita and 
Beek, 1962), The most recent theory has been discussed by 
Astarita et al. (1981). It is assumed that the role of additives 
is that of a catalyst in the reaction rather than a modifier 
of vapor -.liquid equilibrium. Two main types of additives 
have been studied, inorganic and organic; both of which have 
been used commercially. Among the inorganic additives, 
arsenious acid and among the organic additives, ethanolamines 
exhibit large promotional effects. 

The absorption rate promotion by amines has been 
explained in terms of the so called "Shuttle Mechanism" by 
Astarita et al. (1981). In the region near the gas—liquid 
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interface, CO 2 i s assumed to dissolve and undergo a fast 
reaction with free amine resulting in the formation of 
carbamate. This reaction is fast enough to steepen, signifi- 
cant ly, the CC >2 concentration profile, and therefore to 
enhance the mass transfer rate. Figure 2.3 shows a qualitative 
diagram describing the* shuttle mechanism'. The reversion of 
the carbamate ion to bicarbonate ion is known to be a slow 
reaction. Therefore, to react significantly, the carbamate 
ion must diffuse to the bulk of the liquid phase, where it 
reverts to bicarbonate and liberates free amine. The free 
amine then diffuses back to the interface region where it can 
react again with carbon dioxide. The overall reaction for CC^ 
absorption is still eqn.(2,5). But, for rate promotion with 
the help of amino compounds, the reaction is 

C0 2 + RR 1 NH — * RR'NCOOH (2.9) 

(Amine) (Carbamate) 

Reversion of the carbamate to bicarbonate ion takes place 
in the bulk of the liquid as follows: 

RR'NCOOH + OH 1 - — * HCO^ + RR'NH 2 (2.1-0) 

The sequence of reactions given by Eqs, (2.9) and (2.l0)can 
be considered as "homogeneous catqlysis mechanism" or "shuttle 
mechanism" according to the relative rate of reaction (2.10) 
with respect to (2.9). For the "homogeneous catalysis 
mechanism" reaction (2.10) is so much faster than (2.9.) that 
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Fig. 23 Shuttle mechanism for absorption in amine -promoted 
carbonate solution 
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it follows the latter immediately and locally. For the shuttle 
mechanism, however, the reverse is true, i.e. reaction (2.10) 
is much slower than (2.9) and the reaction takes place in the 
bulk of the liquid. At room temperature, the relative rates 
of the two reactions are such that the shuttle mechanism is 
operative. However, at higher temperatures, in the range of 
industrial operation, the rate of react ion ( 2.10) increases 
significantly and the system is better represented by the 
homogeneous catalysis mechanism. This mechanism shoots up the 
overall reaction rate to instantaneous regime. The kinetics of 
rate promotion step will be discussed later (Section 3.5). 



CHAPTER 3 


MATHEMATICAL MODELLING OF HOT CARBONATE ABSORBER 

Since chemical absorption is the basic mechanism of 
acid gas absorption in the Benfield process, the relevant 
theory will be discussed briefly. The phenomenon of gas 
absorption in a liquid with simultaneous chemical reaction 
has been studied extensively over the last three decades 
through wide experimentation as well as theoretical analysis. 
But it was only in the late 60s, that Ramm (1968); Astarita 
(1967) and Danckwerts & Sharma! (T9 published these information 
in a collected, coherent fashion. Since then many works have 
followed (Leder, 1971; Savage et al, , 1980; Astarita and 
Savage, 1980; Astarita et al., 1981 etc.) and the most recent 
analysis and information are available in the treatise by 
Astarita et al. (1983). 

3.1 General Background: 

For the present analysis, we focus attention on a system 
consisting of a chemical solvent used for absorbing one 
component from a gas stream. Let A-^be the component to be 
absorbed which can be present both in gas phase and liquid 
phase. All other components of the gas stream are regarded 
as non-condensable except water vapor. The components present 
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in the liquid stream are considered to be non— volat ile 
except water and the solute A^. Further analysis is based on 
two important terms, molarity of the chemical solvent and the 
degree of saturation or fractional saturation* Molarity, m, 
is the total equivalent concentration of the species that may 
react equimolar ly with component Aj. It is also the actual 
concentration of the reactive species that could be attained 
in the chemical solvent if, in the desorption unit, all 
chemically combined A-^could be stripped from the solvent. The 
degree of saturation or fractional saturation, x, is defined 
as .the ratio of the concentration of the absorbed component to 
the total concentration of the reactive component in both 
reacted or unreacted forms. 

The reactions between A^and the non— volatile 
reactive species in the liquid can be expressed as 



(k = 1,2,...,M) (3.1) 


where M is the total number of independent chemical reactions. 


The stoichiometric coefficient of the jth non-volatile 


components in feth reaction are not independent of each other 
and obey the follcwing equation for preservation of atoms in 
a chemical reaction 

1 + 0(jk = 0 ' k = 1 ' 2 ' — ,M (3 * 2> 



is positive if a component is a reactant, negative if it is 
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a product and zero if it is inert. represents the number 

of moles of A^which are present in chemically combined form 
in component Bj (Astarita and Savage, 1980 ). 

For the case of unpromoted CO^— K^CO^ system there is 
only one overall reaction (M=l ) and there are four species 

to be considered namely A^ = C0 3 (the volatile component), 

2 * 

®1 = OH', = C0 3 and = HCO^ (the non-volatile components). 
The stoichiometric coefficient of Bj and the number of moles 
of A^ chemically combined with B~, denoted by ^ and p.j 


respectively are given as 

follows : 



OH* : 

* 1 * ° 

P =0 


2 ' 

co 3 * 

^2 = 1 

* ^2 = ^ 

( 3 . 3 ) 

HCO^: 

^ = 2 

, P '3 “ 1 



The total concentration of CO 2 in the liquid phase in free 
and combined forms; denoted by a, is given by 

— a, + "b'”' R . b < — a. + b 0 ,+ b_ ( 3 . 4 ) 

1 j J i 2 3 

where a^, b2 and b 3 are the concentration of tine species Aj, 
B2 and B 3 respectively in the liquid phase. 

The liquid phase compass it ion of absorbed component 
described by fractional saturation, x, is given as follows: 



where m, is the molarity, i.e., the total K 2 C0 3 concentra ‘ taon 
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m = [kI /2 = (2b 2 + b 3 )/2 (3.6) 

where [k ] is the potassium ion concentrat ion. 

From (3.4), (3,5) and (3.6) we obtain; 

■ a = -a + m (1+x) (3.7) 


r j t is the rate at which the Kth reaction is taking place 
(for absorption is positive and for desorption it is 
negative) then film theory differential equations for diffusion 
reaction are 



(3.8) 


and 


dy^ T 


(3.9) 


where y is the distance from the interface towards bulk 
liquid and D is the liquid diffusivity, assumed to be same 
for all components. Differentiating equation (3.4) twice with 
respect to y and substituting in equations ( 3.8 ) and (3.9) the 
following equation is obtained; 

— « 

i. M U— V U 

'jk 


2 

d a 


dy 


F[-s < 1 + bh r 1k )] 


( 3 . 10 ) 


In view of equation (3.2), Eq.(3.lO) reduces to 
dy 


(3.11) 


Equation (3.11) implies that the profile of a through the film' 
is linear and therefore 
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d a 
dy 


a o ~ a i 


(3.12) 


where suffixes i and o refer to the interface and bulk liquid 
conditions respectively and 5 is the thickness of the film. 
The physical mass transfer coefficient k° is related to 6 as 

Lj 

k L = "IT (3.13) 

Since BjS are non-volatile, their flux at the interface, when 
all dif f usivities are equal, is zero, i.e., 

db . 


dy 


= 0 


(3.14) 


Combining equations (3.4), (3.12) and (3.14) we arrive at 

(3.15) 


3a 

dy 


a o a i 


If N is the absorption flux of A^at the interface, then by 
definition 

N = - D (3.16). 

dy . 

Equation (3.13), (3.15) and (3.16) when combined yield, 

(3.17) 

A "chemical" mass transfer coefficient is customarily defined 


N = k° ( « ± - a Q ) 


as 


N = k L (a li- a lo> 


(3.18) 


Comparing equations (3.17) and (3.18), we arrive at an important 
concept of enhancement factor, I, defined as, the ratio of the 
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actual rate and the 
same driving force 

I 


rate that would he observed under the 
in the absence of chemical reactions 





(3.19 ) 


In the classical chemical absorption theory, four asymptotic 
regimes have been identified; the kinetic (KR), diffusion al (DR), 
fast reaction (PR) and instantaneous reaction (IR) regimes 
(Astarita, 1967). There are, however, transition regimes like 
kine tic-dif f usional, dif f usional-fast , etc. 

The value of the enhancement factor will depend on the 
particular regime controlling the chemical absorption of C0 2 . 

For reasons to be described, in detail, later in this chapter 
(Section 3.5), we focus attention on the instantaneous reaction 
regime. Physically, instantaneous regime is the case where the 
chemical reactions are so fast that chemical equilibrium 
prevails at all points in the liquid phase and therefore the 
chemical reactions take place at whatever rate the reactants 
are brought together by diffusion. Quite obviously the 
instantaneous regime represents an upper band, since all 
resistance to absorption (or desorption) rate due to chemical 
kinetics has been reduced to zero. Though Savage et al* (1980) 
have given a mathematical expression of the enhancement 
factor, I, for the instantaneous regime, it is still not easy 
to calculate the above quantity as one is required to know 
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the interfacial compositions of A 2/ B 2 , B 3 etc., which in 
practical cases, are very difficult to find. Van Krevelen and 
Hofti jzer (1948) developed an expression for the enhancement 
factor which was basically the Hatta number, defined as 


I 


(_id )V3 ( hl^R ,1/2 

-g-f 2 11 


tanh 


jiia/ 3 r k n C R q/2 

2 ^ ID 

Of 


( 3 . 20 ) 


where p. = viscosity of liquid, Kg/(m)(hr) 
p = density of liquid. Kg/m 

2 

g = acceleration due to gravity, m/(hr) 

3 

k 2I = second-order reaction rate constant, Kmol/(m )(hr) 

3 

C R = average concentration of active component, Kmol/m 


This equationis simple to use and the sophisticated analysis 
of Astarita and Savage (1980) according to the authors themselves, 
have improved very little over the results of the above 
equation (3.20). Estimation of k^ 2 is given in Appendix A, 

3.2 Material and Energy Balances: 

The material and energy balance equations have been 
developed on the basis of the two— film model of mass transfer 
proposed by Whitman (1923). The equation for the flux 
moles of kth component transferred per unit time per unit 
length of the absorber can, therefore, be written as 
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N k = <V a) (p ko- p ki )(s) = < It i k a)(c kl - c k0 )(s) 

” (K gk a> <Pko -P k )(s) (3.21) 

where kg k = gas phase mass transfer coefficient of the kth 
component. Kino 1/ (hr ) (m 2 ) (atm) 
k lk = li( 3 uid phase mass transfer coefficient of the kth 
component, Kmol/(hr ) (m 2 ) (Kmol/m 3 ) 

9 3 

a = wetted area available for mass transfer, m /m 
- partial pressure of component. k in the bulk of 
the gas, atm 

Pki = Partial pressure of component k in the- gas at 
the interface, atm u 

C kQ = concentration of component k in the bulk liquid, 
Kmol/m 3 

= concentration of component k in the liquid 

3 

interface, Kmol/m 

2 

S = area of cross-section of the packed tower, m> 
p* = partial pressure of the kth component in the 
gas phase if it were in equilibrium with the 
bulk liquid phase, atm. 

In spite, of its simplicity, the two-film model, for all 
practical purposes, is a reasonably accurate model. Astarita 
and Savage (1980) have observed that the film theory solution 
not only gives a very reasonable first-order approximation but 
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also coincides with the solution obtained from the more 
sophisticated penetration theory model when all dif f usivities 
are equal. 

The material and energy balance equations are developed 
around a differential height of the packed absorber shown in 
Fig. 3. 1. Envelope III is an elemental volume in the differential 

packed height , dZ , of the absorber, consisting of the gas 
and liquid films denoted by envelopes I and II respectively. 

Nf, are the mass fluxes of the components and dq^ t dq^_ and dq 
are the different heat fluxes at the gas-liquid interface 
separating envelopes I and II. The major assumptions are: 

1. The pressure drop through the bed is negligible. 

2. Steady state conditions prevail. 

Overall mass balance on envelope III: 

dL _ dV 


az - az + rS 

3 

where r = rate of reaction, kmol/(m )(hr) 

L = molar flow rate of liquid, kmol/hr 
V = molar flow rate of gas, kmol/hr 
Z = spatial variable along the height of the 
absorption column, mi 
Overall mass balance on envelope I: 


(3.22) 


dV 

dZ 


= - <K i + V 


(3.23) 
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Interface 
surface = a dZ 



*-,*v x 2>* ,y 2 » T 


Fig. 3-1 Differential section of a packed absorber 
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In the present case, carbon dioxide and water vapour are 
the two components in the g as phase that are absorbed in 
(or desorbed from) the liquid phase* Therefore V*N T consists 
of , the molar flow rate of carbon dioxide per unit time 
per unit length and N 2 , the molar flow rate of water vapour 
per unit time per unit length. 

Component balance on envelope III: 


d(lix k ) d(Vyp 
dZ = ~dZ 


which for the present case yields 


(3.24) 


and 


d( Lx-^ ) 
___ 

d(L>c 2 ) 

dZ 


d(Vy 1 ) 

~dZ 

d(Vy 2 ) 

dZ 


+ rS 


+ rS 


(3.25) 

(3.26) 


where x^. = mole fraction of the kth component in the liquid 
phase in unreacted form 

y, = mole fraction of the kth component in the 
J k 

gas phase . 

Suffixes 1 and 2 denote carbon dioxide and water bapour. 


respectively. 

Component balance on envelope I: 


d(Vy, ) 

= - N k 


(3.27) 


dZ 
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Equation (3*27) gives rise to the following equations for 
the present case, 

d(Vy 1 ) 

dz = *" N l (3*28) 

d(Vy 2 ) 

“dZ = - N 2 <3.29) 

Combination of the equations (3.23) and (3.28) gives 

dy l _ N 2 y l ** 

dZ = V (3.30) 


Similarly equations (3.23) and (3.29) yield 

a y 2 _ “1*2 - N 2 (1 -y2> 

dZ V 


(3.31) 


Again, combination of equations (3.22), 

gives dx-^ (rs— ) (l-oc^ ) + X ]_ N 2 

___ = _ 

Similarly, from equations (3.22), (3.23) 

arrives at the following equation 

dx 2 (rs~N 2 ) (1-k 2 ) + x 2 N l 
dZ“ = L 


(3.23), (3.25) and. (3.28) 

(3.32) 

, (3.26) and (3,29), one 


(3.33) 


The energy balance of the packed absorber is done on the 
consideration that it operates under adiabatic condition 
* Heat balance on envelope II: 

By making energy balance for the liquid side around 
the differential tower height and taking help of equation (3.22) 
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one arrives at 

*L = _ (H l+ N 2~r S)t ' _i_ to i ' 1 

aZ L " L C PL 32 " ^pL SZ ‘ ^pL 

(3.34) 

where dq^ = rate of heat carried by the absorbed components 

from gas phase to liquid phase over the differential 
height , dZ. 



= (N 1 C pl + N 2 C p2 ) dZ T, Kcal/hr (3.35) 

where C p ^ and C p2 are the average specific heats of carbon 
dioxide and water respectively in the gas phase and C pL is the 
average specific heat of the liquid. The values of C ^ and C p2 
are 8.1072 and 9,9023 Kcal/(Kmol) (K) , respectively. 

dq t = rate of heat transferred from gas phase to liquid 
phase in the differential section of the tower due 
to temperature difference over the differential 


height, dZ 

= hadZ.S (T— t ) , Kcal/hr (3.36) 

where h = coefficient of heat transfer in the gas-liquid 
interface, Kcal/(m^ ) (hr ) (k) 

T = temperature of the gas phase, K 
t = temperature of the liquid phase, K 
dq = rate of heat generated due to absorption in the 
differential height dZ of the tower 

(3.37) 


(rS + N 2 4H 2 )dZ, Kcal/hr 
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where = heat generated due to react 


ion of carbon 


dioxide with hot potash, Kcsl/Kmol 
^ H 2 = hea t of vapourisation of water, Kcal/Kmol 

Incorporating the expressions for dq , dq. and dq (equations 

g 

(3.35), (3.36) and (3.37)) in equat ion( 3. 34 ) one arrives at 


(N l +N 2~ rS) (N l C ol +N 2 C -n? )T 


haS ( T«»t ) 
LC j 


rS 4 h 1 +n 2 

~ LC _ 


(3.38) 


Heat balance on envelope I: 

The gas phase differential heat balance equation is 


given as. 


W _ J^(%) _ _L_ & 


VC v dZ 


VC dz^ 


(3.39) 


Pg Pg- 

Combining equations (3.35), (3,36) and (3.39) one gets 
dT (N 1+ N 2 )T < n i c pl + N 2 C p2 )T haS(T-t) 


(3,40) 


The change in compositions of and KHCO^ are related to 

the change in C0 2 composition according to the stoichiometry 
of the overall reaction (2.2). For the special case of 
instantaneous reaction regime, where chemical reaction takes 
place at the same rate at which the' reactants are brought 
together by diffusion, N-^rS and equations (3.22), (3.32), (3.33) 
and (3.38) will take the forms as follows 


dZ 2 


(3.41 ) 
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and 


x 1 N 2 

dZ L 


dx 2 _ Ni.N 2 (1.x 2 ) 
dZ L 

dt = V (N lSl +N 2 C p2 )T 

dZ L LC _ 

PL 


(3.42) 


(3.43) 


haS ( T-t ) N 1 ^ H l +N 2 h 2 


LC 


pL 


LC 


PL. 


(3.44) 

Equations (3.23), (3.30), (3.31), (3.40), (3.41), (3.42), (3.43) 
and (3.44) form the model for the absorber which must be 
solved together to obtain the unknown parameters. The solution 
procedure is discussed in the following section. 


3.3 Computational Method: 

In a countercurrent absorber , the rich feed gas enters 
at the bottom and absorbing lean liquid enters at the top of 
the column and are usually completely specified. However, 
the lean gas leaves from the top and rich liquid leaves' at the 
bottom, both of which are, at best, partially specified. 
Therefore, the state of neither .end of the column is known 
completely which leads to a two-point boundary value problem 
when integrating heat and mass balance equations through the 
tower. Equations (3.23), (3.30), (3.31), (3.40), (3.41), (3.42), 
(3.43) and (3.44) give rise ttxa highly coupled system of 
equations. For the present case, the packed height was divided 
into three sections as shown. 
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Z - W- t 
iz 


Part of the lean liquid, L, enters at the top and remaining 
liquid, L >2 enters after first section of the packing. The 
lower section is again divided in two parts II & III, to 
enable redistribution of the liquid. In the present scheme, 
each section has been divided in 20 equal integration steps 
giving a total of 60 segments for the total column. This 
step size is found to be satisfactory since further reduction 
in size did not significantly change the results. Over each 


integration step, the thermophysical properties were assumed 
to remain constant. These properties were updated at the end 
of each integration step by procedures to be discussed in the 
next section. 

The boundary conditions known at the two ends of the 
absorber are the 

L( Z=Hj*) = L fc 

x 1 (Z=H t l)= x lt (3.45,) 

x 2^ Z=H t x 2t 
t (Z=H t )= t t 
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V (z=0) 

= v b 

y 1 (z=o) 

11 

y 2 (z=o) 

11 

■< 

T (z=o) 

= T b 


where suffix t denotes t ewer top and b, tower bottom and H^is 
the total packed height of the tower. 

Several methods are available for solving two point 
boundary value problems. For the present work, "Shooting 
method" has been used (Fox, 1957). In the shooting method, 
values for the unspecified conditions at the initial point of 
the interval ("missing initial conditions") are assumed, and 
the differential equations are numerically integrated to the 
terminal point by applying essentially an initial value algorithm 
("Shooting" at the target terminal points). If the computed 
terminal values satisfy the specified terminal conditions, the 
problem has been solved. If they do not,, the differences 
between the computed and specified terminal conditions (the 
"miss distances") are used to adjust the assumed initial 
conditions. If the differential equations and boundary 
conditions are linear, the adjustments need only be made once, 
but if the differential equations are nonlinear, the adjustment 
of the missing initial conditions is an iterative process. 
(Roberts and Shipman, 1972; Fox, 1957; Keller, 1968). In the 
present case, the differential equations are 
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coupled. The integration is started from the bottom of the 
tower with two sets of assumed values for the unknown variables 
at the tower bottom, i.e., L^, x^, x 2 - b and t^ # as well as 
related concentrations of other components in the liquid 

x 3 b' x 4b / X 5b / x 6b ^ or K 2 C0 3 / khc ° 3 ' amine and V 0 0,_ respectively. 
The integration is performed using fourth— order Runge—Kutta 
method. The known boundary conditions at the tower top, i.e. 

^t • x lt / x 2t an< ^ ^t ' are compared with those computed for 

the two assumed sets. The difference between the known and 
computed values for the boundary conditions at the top form 
the difference functioning. 

Figure 3. 2 shows how the shooting method works. For 


the liquid flew rate at the bottom, two values and L^ 

are assumed. The 'corresponding calculated values ca 2 _ c and 

L, _ 1 are used to calculate the difference functions 

t2^calc. 

L t - L tl,calc. and L t - L t2,calc. < where L t 13 the observed 
liquid flow rate at the top) which are plotted as A and B. 

Line AB is extended to intersect x-axis at C which provides 


the next guess for 1^ say L^. For the next iteration, Lj^ and 
are the two values used to get the new guess in the same 
fashion. The remaining 7 variables, x ib' - x 2b / x 3b' x 4b' x 5b' 
x 6b an< ^ ^b 50:6 '* :rea ^ e< 3 an ^he Same manner. 

A computer program has been written in FORTRAN language 
for DEC 1090 system using above procedure. A brief User's 
Mannual is included in Appendix b. 
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Before we can actually integrate equations (3.23) to 
(3.44) we need to obtain the thermodynamic properties such 
as partial pressure, Henry's law constant etc. and transport 
properties such as viscosity, thermal conductivity and so on. 

These are discussed in the next section. 

3.4 Thermodynamic and Transport Properties : 

Evaluation of the requisite thermodynamic and transport 
properties constitutes a major part of any mathematical 
modelling effort. Use of appropriate property estimation 
equations is an integral part of good design program. The 
original data on vapour-liquid equilibrium for CC^-hot carbonate 
system were published by Benson et al« (1954; 1956), Tosh et al. 

(19 59, 1960). The other physical properties of the system 
such as density, viscosity, specific heat etc. of the hot 
carbonate solution were published by Danckwerts & Sharina(1966 ) ;Bocard| 
and Maryland (1962), Ratcliff and Holdcroft (1963), Kohl and 
Risenfield ( 1979 ) and others. More recently, Kent and Eisenberg 
(1975), Savage et al. (1980), and Astarita et al. (1983) have 
published latest data on thermodynamic and transport parameters. 

For the gas phase, which is at lew reduced pressure and at high 
reduced temperature, ideal thermodynamics can be used since 
the compressibility factor is close to unity. In the present 
case, the criterion is reasonably satisfied and herce ideal 
thermodynamics have been used to estimate various properties. 
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ihe liquid phase properties have been estimated using available 
correlations or curves as described below, 

3,4,1 Vapour— Liquid Equilibrium (VLE ) : 

For the system of C0 2 in. K 2 C0 3 /KHC0 3 solution, a single 
parameter can be defined which completely characterizes the 
chemical composition of the liquid phase. The parameter chosen 
for this case is the fractional saturation, x. For an absorption 
system, x increases along the liquid path in the absorber. If 
p* represents the partial pressure of the absorbed component 
which equilibriates the liquid phase composition; the VLE 
relationship becomes (Astarita et al, 1983) 

p* = Am < ^F(x / t) (3.47) 


A 1 signifies a parameter into which all the activity coefficients 
are lumped. In spite of this gross simplification, the above 
equation is usually quite effective in predicting the correct 
form of VLE relationship. The value of the exponent q in 
equation( 3.47 ) depends on the stoichiometry of the main reaction 
involved and F(x,t) is some function of x and t. Actual data 
on C0 2 ^K 2 C0 3 system by Tosh et al.(1959), Bocard and Maryland 
(1962) and also by Kent and Eisenberg (1975) are well correlated 
by the following equation suggested by Astarita et al. (1983). 


p* = 1.95 x 10 9 m? ,4 ( 


x 2 . , -8160b 

) exp ( —— ) 


1-x 


t 


(3.48) 
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wnere pj, the equilibrium partial pressure of C 0 2 , is in 
atmospheres, m is expressed in Wnol/m 3 and t in K. The 
constant A 1 in Eq.(3.47) taken here as 1.95 x 10 9 , can vary in 
case of very high molarity and especially when the effect of 
promoter is strong. 

For the water vapour a different equation was proposed 
by Tosh et al. (1959). On the basis of the assumption that the 
water vapour obeys ideal gas law. The following Clausius- 
Clapeyron form of equation was put forward by them and was 
found to correlate their data very well. 

log p* = - T Tj 3 03xR ) ( £- : ~ 2. 45x1 O * 3 )— 1 .1672+C 0,49) 

where p* is the equilibrium partial pressure of water vapour 
in atmospheres? L is the heat of vapourization of water, and C 
is a constant for a constant molarity or weight percent of K 2 C 0 3 
solution. The' heat of vapourization was taken to be 9,800 Kcais/ 
K mole which is close to the average value between 70°C and 
130°C obtained from the steam table. The value of C decreases 
as conversion of K 2 C 0 3 to KHC 0 3 increases and the following 
expression has been fitted to the data of Tosh et al, by 
nonlinear regression 

C ~ a c + k G / x c + V x c 

where x Q is the conversion of K 2 C 0 3 defined as 

0.691 (% KHCO 3 ) 

x c = % K 2 C0 3 + 0.691 (% KHC0 3 


(3.51 ) 
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and a G , b c and are regression coefficients. 

The total equivalent weight percent K 2 CC> 3 is defined as 


Equivalent K 2 C0 3 % = 


%K 2 C0 3 + 0.691 (% KHCC> 3 ) 

'?k 2 C0 3 + 0.781 (%ICHC0 3 ')+9il 0 X 100 


(3.52) 

For a 30% equivalent K 2 C0 3/ a c , b c , and were found to be 

a c = 1*2014 
b Q = 0.23757 
C G = -0.0537 

3.4.2 Henry's Law Constant: 

The experimental measurement of the value of Henry's 
law constant, H, is quite straight forward for a physical 
solvent; whereas, for a chemical solvent the situation is 
considerably more complex. Moreover, for an electrolyte solution, 
the ionic interactions tend to affect its value in a complex way. 

If H° represent the value of the Henry's law constant for pure 
solvent of a chemically reactive solution (e.g., water, in the 
case of aqueous solutions), the value of H, for the reactive 
solution, will usually be of the same order of magnitude as H° 
but appreciably different from it. This is due to the fact j 

that appreciable concentration of the reactive component, 
usually an electrolyte, increases the "salting out" effect. 

Van Krevelen and Hoftijzer (1948) and later Danckwertst&’ Sharma(l9 66: 
put forward, on theoretical basis, a correlation as follows: 
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l0g ^ = K i x i (3.53) 

where 1^ 

where Z. 

C. 

i 

and K ± 

For the CO 2 — K 2 C0 3 system, equation (3.53) applieswith 
3 

= 0,062 m /Kmol. For this system, the ionic strength can 
be related to molarity as 

1^ = m (3-dc) (3.54) 

and an alternative expression for H in terms of molarity can 
be written as (Astarita, et al., 1983)? 

log — = K m. (3.55) 

H° m 

3 

where = 0.025 m /Kmol for C0 2** K 2 C0 3 system. 

The above equation correlates well, the data given by Kent and 
Eisenberg (1975) and Astarita et al. and has teen used in the 
present model. 

The Arrhenius plot for inverse Henry* s law constant of 
C0 9 in water, i.e. , [h 0 }* 4-1 given by Astarita et al.(1983) has 
been correlated, for the present work, by linear regression 


- the ionic strength of solution, Kgion/m 3 

= i 21 Z? C. 

2 ll 

= charge of each ionic species 

= concentration of each ionic species, Kg~ion/m 3 
= a constant that depends only on the nature of 
the ions present in solution, m 3 /Kmol. 


analysis as 
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log [h ] _ ajj + bj^/t (3.56) 

where a^ and bj^ are regression coefficients and have the 
values —4.3856 and 867,4932, respectively. 

3,4.3 Mass Transfer Coefficients: 

Several correlations are available for estimation of 
gas— side mass transfer coefficient .Van' Krevelen and Hoftijzer 
(1948), Danckwerts (1970), and Uchida et al,(l972) have reported 
a number of such correlations^ Astarita et al,. (1983) have 
recommended the following correlation proposed by Onda et al. 
(19 68), which has been used in the present model. 



where the subscript g refers to properties in the gas ;phase, 

G is the mass flow rate of inert gas in Kg/Cm )(hr), d is the 

packing size, and C* is a dimensionless constant which depends 

on the size of packing. For the present case," where 1.5 

Intalax saddles have been used, C' is taken as 5.28. R is the 

universal gas constant in (m ) ( atm)/(Kmol ) (K) , ^*g ' Jg anG 

D are the viscosity, density and diffusivity of the kth 
gk 

component respectively in the gas phase. i n this case d p =O.Oll5. 

The liquid side mass transfer coefficient is obtained 
from the analysis of Van Krevelen and Hoftijzer (1948) as 


follows : 
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*^> V3 


sJL 


D 


0.015 (-i-)V 3 _a£ 

all' 'HfE' ,,2 




tanh l^^ihfS-yk 


g r 


D 

(3.58) 


"where L Xic[uic3. flow rete^ Kief,/ (m )(sec) 

g = acceleration due to gravity, m/sec^ 

f = liquid density. Kg/m 3 

l-l = liquid viscosity, kg/(m)(sec) 

2 

D = liquid diffusivity, m/sec 

3 

kjj= reaction rate constant, kmol/(m )(sec) 

= average concentration of active component, 
kmol/m 3 

Experimentally it was found that the above semiempirical 
equation is a satisfactory approximation for the correlation 
of the data. 

The overall mass transfer coefficient for the absorbed 
component CC> 2 consists of the gas and liquid side mass transfer 
coefficients and is expressed as 


K 


•gi 


_i_ + H. 

k gl hi 


(3.59) 


where H is obtained from equation (3.55). 

It may be assumed that there is no liquid side resistance xor 
£he mass transfer of the solvent water vapour. Therefore, the 
overall mass transfer coefficient for the solvent water vapour 


IS 
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K - = k - 
g2 g2 


(3.60) 


3,4.4 Heat Transfer? Coefficient; 

Due to lack of experimental data on gas-liquid heat 
transfer coefficient in packed column, Uchida et ai.(1972) have 
recommended that the analogy between the heat and mass transfers 
can be used to calculate heat transfer coefficient, h_, in the 
packed column. Thus 

Jfj — jjj (3.61 ) 

where j H = j factor for heat transfer 
and j-Q = j factor for mass transfer 


By definition. 


,• _ h ( ^g Efl ) 2 / 3 

J H C f v A 
pg ^g 


(3.62) 


and 


i _ ! ( jjg _ ) 2 / 3 

4d " v Tg D gk 


(2.63) 


where suffix k denotes the kth component transferred from 
gas to liquid. 

Combining (3.61), (3.62) and (3.63) we get, 

2 £ 

h =k ( A g pg— )V 3 0.64) 

wheire C = specific heat of gas mixture, Kcal/(Kg)(K) 

3?Q 

J\^g = thermal conductivity of gas mixture, Kcal/(hr ) (m ) (K) 
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Equation (3.64) has been used in the present model to 
compute component heat transfer coefficients. Overall heat 
transfer coefficient, h, for the gas mixture has been obtained 
by taking the molar average of canponent coefficients. The 
experiments performed by Heertjes and Ringeps (1956) and 
Linch and Wilke (1955) show the validity of this relation. 

3.4.5 Wetted Interfacial Area: 

Usually, interfacial area reported for various sizes of 
the dry packing is different from the gas-liquid interfacial 
area of the wetted packing, that is actually available for 
mass transfer. Van Krevelen and Hoftijzer (1948), Danckwerts 
(1970), Puranil and Vogelpohl (1974) and Laurent and Charpentier 
(1974) have discussed the estimation of wetted interfacial area. 
For the sake of simplicity we use the correlation proposed by 
Van Krevelen and Hoftijzer given as follows: 



2 3 

where a^ = dry packing surface area, m /m 

2/3 

a = wetted interfacial area, m /m 

2 

L = liquid mass flow rate, Kg/Osec )(m ) 

3 

f = density of liquid. Kg/m 
C = 5000.0 

0L 

The same interfacial .area has been used for heat transfer 


3ilS O* 
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3.4,6 Density of Liquid and Gas: 

Bocard and Maryland (1962) and Kohl and Risenfield (1979) 
have reported specific gravity data of aqueous solutions for 
different equivalent weight percent K 2 CC> 3 at different 
temperatures. Following relation has been obtained by linear 
regression of the data of Bocard and Maryland both with respect 
to temperature and concentration, which is valid for temperature 
range of 70 to 130°C and 20 to 40 weight percent total equivalent 

K 2 c °3 . 

s = a + b 0 w <- c- t„ (3.66) 

So O V' 


where s = specific gravity of K 2 C0 3 so l ut; j- ori 

w = total equivalent weight percent K^CO^ 
tg = liquid temperature, °C 

a g/ b s and c g are regression constants with the following 
values 

ag = 1.0679 
b s = 0,01 
Cg = 9.7xlO~ 4 

The density of the gas mixture is calculated from the 


ideal gas law as 




P M 


-SL 


RT 


where P = density of the gas mixture. Kg/m 
J g 

P = total tower pressure, atm 
M = molecular weight of the gas mixture 

g 


(3.67) 
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R - universal gas constant = 0.08205 m 3 atm/(K) (Kmol) 

T = gas temperature, K 

3.4.7 Viscosity of Liquid and Gas : 

Data on the viscosity of aqueous K 2 C0 3 - KHC0 3 solutions 
has been given by Bocard and Maryland (1962). The variation of 
viscosity has been reported both with respect to temperature 
and weight percent K 3 C0 3 total equivalent. For the present 
case of 30 weight percent K 2 C0 3 , the following relationship 
has been developed by non-linear regression and used in the 
model 

jj, = a* + b't p + c't 3 + d't 3 (3.68) 

where p, = viscosity of liquid in Kg/(m)(hr) 
tp = liquid temperature, °F 

and a'jb'/C* and d’ are the regression constants with the 
following values : 
a* = 12.9275 
b* = -0.1114 
c* = 4.1184 x 10- 4 
d' = -5.796 x 10~ 7 

The viscosities of the individual gas components at 
atmospheric pressure f averaged over the temperature range of 
70 ~ 135°C have been obtained from Perry (1973). The mixture 
viscosity is obtained by using the correlation of Wilke (Reid 
et al #/ 1977) given as follows: 
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n 


Y i ^ qi 
n 


(3.69) 


1=1 y 0 

J ij 
j=l 


where ti g = viscosity of the gas mixture 

^gi an?l ^ gj are the viscosities of ith and jth 
components respectively in the gas phase 
y^ and y^ are the mole fractions of ith and jth components 


an d 


hj - 


respectively in the gas phase 

_i+( 2 

8(1 + M./M .) * 

l J 


(3.70) 


i and j refer to ith and jth components respectively. 
Since at high reduced temperatures or low reduced pressures 
the viscosity of a gas does not vary significantly from its * 
low pressure value, there is no need to incorporate pressure 
correction for the mixture viscosity given by equation (3.69). 


3.4.8 Dif fusivities of Liquid and Gas: 

Danckwerts and Sharma (1966) and Ratcliff and Holdcroft 
(1963) provided early correlations of diffusivity of dissolved 
carbon dioxide in solvents at different temperatures. Savage 
et al. (1980 ) have given more recent data on diffusivity 
of values for C0 2 ~K 2 C0 3 system. For each of computer calculations, 
the following correlation has been developed by means of linear 
regression analysis of their data 

logD = a^ - bj/t - 0.4437 (3.71) 
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where D liquid diffusivity (assumed equal for all components ), 
m 2 /hr 

a D and b D are regression constants with following values 
a D = -3.0188 

b D = 586.9727 

The binary gas dif fusivities are obtained using the following 
correlation proposed by Sherwood and Pigford (1952) 

2/2 

D kj = (2.6544x10 ) p(v l/3 +v l/3 ) 2 / + / (3 * 72) 

^ j 

2 

where = binary diffusivity in gas phase, m /hr 


V k' V j = m °l ecu l ar volumes of k and j 
The mixture diffusivity has been calculated using the simplified 


correlation suggested by Reid et al. (1977) as follows: 

n 

— ,i 


V = < 


> 


D, 


•) 


(3.73) 


j/k 


kj 


where D , is the diffusivity of the kth component in the gas 
2 

mixture, m /hr. 


3.4.9 Specific Heait of Liquid and Gas: 

Data on the specific heat of K 2 C0 3 solution were 
published by Kohl and Risenfield (1979). The accuracy of these 
data has been reported to be + 0.0024 Kcal/Kg K. These data 
have been well correlated by the following equation and used 


in the present model. 
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C pL “ a " + b " w + c’V 2 + d"w 3 (3.74) 

wnere C pL = specific heat of liquid, Kcal/(Kg)(K) 

a- / b , c and d" are regression coefficients with values, 
a" = 0.9998 


b" = -0.00914 

C" = -0.1063 x 1Q~ 3 

d" = 0.3058 x 10~ 5 


The value of the specific heat of liquid in Kcal/(Kmol)(K) 
i.e. C pL , is obtained by multiplying by the molecular weight 
of the liquid solution. Average specific heats of the individual 
gas components for the temperature range 70 - 135°C have been 
obtained from the nomograph given in Perry (1973). 

The specific heat of the gas mixture has been, calculated 
assuming the mixture to be an ideal one by 


pg 


n 

= 21 

i=l 


y. C . 
J x pgi 


(3.75) 


where C ^ = average specific heat of ith component in the 
gas mixture, Kcal/Kmol/K 

3.4.10 Thermal Conductivity of the Gas Mixture: 

The average thermal, conductivity of the individual gas 
components for the temperature range 70°C — 135°C have been 
obtained from ttie nomograph of Chang (Reid et al., 1977). The 
thermal conductivity of the gas mixture is obtained by using 
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the TJassil jewa equation (Reid et al„ 1977) given belt™ 


w here \ 


gi 





i=i 


Y J. Agi 

m 

V y. a. . 
j=l 


(3.76) 


thermal conductivity of pure component i in the 
gas mixture 

is found by theoretical analysis to be same as 
jZhj equation (3.70) 


3.4.11 Heat of Absorption 

Absorption of CO 2 in K^CO^ s °l ut i°n I s acconpanied by 
two heat effects, the heat of solution and heat of reaction., 
The net heat of absorption has been correlated from the data 
of Bocard and Maryland (1962) by non-linear regression. 
Following relation has been obtained. 

= a h + bjW + c h w 2 + d^w 3 + e h w 4 (3.77) 


where = net heat of absorption of 

hot potash solution, Kcal/Kmol 
a h , b h , c^, 6 ^ and e h are the regression 


the following values 
a h = 6510.8058 
b h = 19.4580 
c h = 1.2727931 
= -0.0313384 - 
e, = 3.955 X'10 4 


carbon dioxide in 


constants with 
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3*5 Effect 


of Rate Promotion on VIE and Kinetics: 


In Section 2.3 it was discussed how the rate of CO 
absorptionis promoted by addition of amines to the K 2 C0 3 
solution. These amines affect the vapour— liquid equilibrium 
and the kinetics of the overall reaction in a complex way. 

The amine can act as a base or it can form carbamates for 
which the chemistry is even more complex. For a base promoted 
solution, where carbamate formation is either neglected or 
considered to be an unstable process, the protonation constant 
of the base B is given as 


= 


[B] [H-*] 


(3.78) 


[ BH + ] 

where [b] = concentration of free base, Kmol/m 3 

[h 3 = concentrationof hydrogen ion, K ion/m 3 
[BH*] = concentrationof protonated base units, Krnol/m 3 
The amine is a weaker base than the carbonate and hence its 
protonation constant Kp is larger than the second dissociation 


constant, K of carbonic acid 


K, 


[ H + ] [CO"] 
[hco»] 


(3.79) 


where [co"l = concentration of carbonate ion, K-ion/m 
3 

[HCO^ ] = concentrationof bicarbonate ion, K-ion/m 
Therefore, the following reaction 
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B + HCO* BH + + CO^ ( 3 .80 ) 

is oisplaced to the left. Hence as long as x does nob approach 
the value of unity, the base will remain essentially in its - 
undissociated from B (Astarita et al., 1983). Reaction (2.5) 
is still the main reaction, and the value of will essentially 
be same as in an unpromoted solution given by equation (3,48). 

At values of x very close to unity, the [HCO^] /[C0^] ratio 
will become quite large and BH + will begin to form, i.e., the 
base will begin to exhibit its ability to react with carbon 
dioxide according to the following reaction 

B + C0 2 + H 2 0 ^±:BH + + HCC> 2 (3.81) 

Values of x in excess of unity can be achieved in this case. 
Astarita et al. (1983) have proposed the following rough rule 

of thumb. 

i) When x <£. 1, the main reaction is (2.5) and the 

base (amine in present case) acts simply as a catalyst. 
The VIE relationship is unchanged and given by 

eqn.(3.48). . 

ii) When x 1, the main react ionis (3.81) and the base _ 
acts as a chemical sink for CC> 2 and the VIE relation- 
ship has to be changed as follows: 

* _ - (x-1) '(l+xl (3.82) 

P “ K„ C + 1-X 

C 1 
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where me is the total amount of amine added to the solution 
and K c _^ is the first ionisation constant for carbonic acid 
and is given by 


*0 


1 


[H + ] [ HCCd ] 

[co 2 ] 


(3.83) 


I£n the present case, x varied from about 0.2 to 0.67 along 


the liquid path of the tower and hence we can safely neglect 


equation (3.81) and consider case i) to be valid. Therefore, 

VLE relationship given by equation (3.48) holds good. 

To see how the kinetics is affected by. the presence of 

amine, we consider the shuttle mechanism once again (discussed 

earlier in Section 2.3)., The shuttle mechanism involves the 

following two chemical steps, 

C0„ + Promoter -> intermediate (3.84) 

Z 

Intermediate + OH' > HCO^ + Promoter (3.85) 

Since these two reactions take place at the interface and in the 
bulk respectively they are, therefore, separated by the diffusion 
processes. In absence of diffusion processes, the sequence 
(3.84) and (3,85) is a classical sequence of homogeneous catalysis. 
Therefore the amine promoter that acts via a shuttle mechanism 
at low temperature may well act as a homogeneous catalyst at 
higher temperature, where reaction (3.85) may be fast enough 

to take place at the interface region. 

If the amine promoter is regarded as a homogeneous 

catalyst as suggested in case i), its largest possible effect 
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is to speed up the rate of reaction (2.5) to the point where 
mass transfer takes place in the instantaneous reaction regime. 

As the ratio of diffusion time to reaction time is much greater 
than unity (about 10^), far the unpromoted solution in the 
present case, the reaction is fast enough to result in signifi- 
cant rate enhancement; and, on being promoted, the rate promotion 
factor E is most likely to reach its maximum value, defined as 




E 


max 


(3.86) 


unpromoted 

where I ^ is the enhancement factor in'the instantaneous regime 

and I unpromct ed is that f or the unpromoted solution. 

Hence the effective rate enhancement factor for the present 

case is 


"’'promoted = max X unpromoted “ 


(3.87) 


The value of 1^ has already been given by equation (3. 20) 



CHAPTER 4 


RESULTS AND DISCUSSION 

4*1 Calcu lated Results and Model Validations 

The model presented in the previous chapter has been 
tested for carbon dioxide absorber in an existing ammonia 
plait . The plant uses hot potassium carbonate solution with 
diethanolamine (DEA) promoter for absorption of carbon 
dioxide. The packed absorber has three sections. The lean 
solution from the regenerator is split into two parts. One 
part enters the absorption tower at the top after being 
cooled , whereas, the other part enters after the first section 
of the packing (see figure in Section 3.3). Table 4.1 shows 
the design parameters for the absorption tower. Table 4.2 
shows the value of the input and output variables such as 
flow rates of liquid and gas streams, the composition of each 
stream, temperature etc. Also shown in Table 4.2 are the 
calculated values for the same outlet variables for comparison 
Except for the mole fraction of water vapour in the 
gas phase, all the calculated liquid and gas compositions 
match well with the plant values very well, the maximum 
deviation being 7 percent for the carbon dioxide. Even this 
7 percent deviation may be due to measurement error since 
there is a higher uncertainty in measurement for small 
concentration range. If we calculate the percent co 2 remov 
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Table 4.1: Design Parameters for the Absorber 


Height of Packing : 

i ) 1st sect ion : 

ii) 2nd section : 

iii) 3rd section : 

Dia. of packed bed : 

Packing size and material: 


6*0 m 
5.6 m 
5.6 m 

1.4 m 

37 mm. (1.5") Ceramic 

Intalox saddles 
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Table 4.2: Comparison of Calculated Results 
with the Observed Plant Data 


Par ame ters 


Inlet 


Outlet 


Observed Calculated ^Deviation 


Liquid flow rate, 
Km o 1/hr 


Liquid temperature, 

K 


2629.8357 

343.0 


5308.0 5351.690 *.0.82 


Liquid composition : 
i) Mole fraction of 


K 2 C0 3 


ii) Mole fraction of 
KHC0_ 


0.04013 

0.02015 


iii ) Mole fraction of 
water 


iv ) Mole fraction of 
DEA 


0.9320 
7.193x10' 


-3 


v) Mole fraction of 


V 2°5 


Gas flow rate, 
Kmol/hr 


5.554x10' 

762.877 


-4 


Gas temperature , K 408.0 


Gas composition: 

i) Mole fraction of 


CO, 


ii) Mole fraction of 
CH , 


0.1659 

0.00184 


iii) Mole fraction of 
CO 


0.00263 


iv) Mole fraction of 


H. 


v) Mole fraction of 

n 2 

vi) Mole fraction of 
AT 

vii) Mole fraction of 


0.5306 

0.1733 


0.0021 

0.1237 


392.0 

390.0 

0.50 

0.0158 

0.0159 

-0.63 

0.0663 

0.0669 

-0.90 

0.9198 

0.9095 

1.12 

7 . 0053x1 0“ 3 7 .070x1 0~ 3 

-0.92 

5, 39 9x1 0~ 4 5. 449x1 0* 4 

-0.93 

547.090 

544.69 0 

0.43 


382.2 

- 

0.001 

0.00107 

—7 .0 

0.00257 

0*00258 

*0.39 

0.00366 

0.00368 

—0. 55 

0.7388 

0.7431 

—0 . 58 

0.2413 

0.2427 

*0.58" 

0.0029 

0.00292 

*0.69 

0.0099 

0.0058 

41.40 
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then the observed value is 99 4 ^ „ 

Xb yy.4 Which is to be compared with 

the calculated value of 99.36. The calculated value of mole 
fraction of water vapour in' the gas was 41.4%. lower than the 
plant data. While the actual reason for this discrepancy is 
unknown, it may be because of the vapour-liquid equilibrium 
data (Tosh Gt al . , 1959) which are old and may be in error. 

Again if wo consider the percent removal of water vapour in 
the column then the observed value is 92% which must be compared 
with the calculated value of 95.3% and now the match looks 
more reasonable . The flow rates and temperatures of both gas 
and liquid streams match well. Since most of the data match 
quite closely with the calculated values, it can be concluded 
that the model is satisfactory. 

To establish the generality of the model developed in 
the present work, it is necessary to validate it against several 
nets of plant delta preferably in different regions of temperature, 
flow rates, etc. However , getting such actual data is not 
easy and in the present case, only one set of data could be 
obtained. 

4.2 Calculated Profiles: 

The temperature, flow and concentration profiles have 
been calculated using the model along the packed tewer hezght 
from the bottom to the top. A typical set of profiles are 
shown in Figures 4.1 and 4.2 for the data set given in 





Concentration, mole fraction 



Cone of H2O in liquid, mole fraction 
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Tables 4.1 and 4.2. The step change in liquid flow rate and 
liquid temperature seen in Figure 4.1 is because of the 
addition of split stream of lean absorbent at the end of first 
section of packing. While the gas flow rate is not affected 
by the liquid addition, the gas temperature shows a slight 
increase. Since the ted height is measured from bottom up, 
the addition of liquid shows a step decrease when the x coordi*- 
nate changes frcm 0 to 1. Gas flow rate decreases along the 
height because of absorption of CO^ and water vapour. The gas 
temperature decreases as it moves up the column since heat is 
being transferred to liquid because of absorption and temperature 
difference. For the same reasons, liquid temperature increases 
as it moves down the column. 

The concentration profiles fcr various components are 
shown in Fig. 4.2. The concentrations of CC> 2 and water vapour 
in gas phase decrease as the gas moves up the column. Since 
CO 2 reacts with K 2 C0 3 there is no free carbon dioxide in the 
liquid phase. Concentration of ^ 2 C0 3 in fche solution decreases 
as it moves down the column since it reacts with the dissolved 
CO2. As the reaction product is KHCO^ , its concentration 
increases as the liquid travels towards the bottom. 

4.3 Sensitivity Analysis; 

From a plant operator's point of view, it is important 

the absorption process with respect 


to know how sensitive is 
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to certain operating variables, in the present case where 
the objective is to remove CC^ frcm the synthesis gas, one 
would like to find hew the exit stream CO^ concentration 
would vary if the operating pressure of the tower was to 
change from its design value. Figure 4.3 shews this result. 

A decrease of 5 atms in total pressure could significantly 
lower the CO^ removal efficiency. Figure 4.4 shows how exit 
CO 2 concentration would change if the fraction of lean 
absorbent fed at the top of the column is varied from its 
design value. This is a less sensitive parameter compared 
to the total pressure, but, an increase in this stream will 
result in a higher purity of synthesis gas. Changing the 
temperature of cold lean stream fed at the top of the tewer 
seems to have very little effect on CO^ removal (Figure 4.5, 
curve A). This is because in the low temperature range, VIE 
is not significantly affected by temperature variation. Also 
the amine promotion is negligible at low temperatures. However 
an increase in temperature of the split stream lean absorbent 
adversely affects the gas purity as seen from curve B in the 
same figure* 








Temperature ,°C 


Exit C0 2 concentration Vs. temperature of lean streams 




CHAPTER 5 


CONCmsiONS AND REC QMMENEAT IONS 

The model of Cc >2 absorber, developed in the present 
study, adequately accounts for the effect of presence of 
amines in the hot carbonate solution. While amine seems to 
play the role of a catalyst and pushes the C0 2 ~K 2 CC> 3 reaction 
to the instantaneous regime, the promotional effect in modelling 
the VLE curve is negligible. A close agreement between the 
calculated values and experimental results for one set of 
plant data justifies the assumptions made.- However, in 
absence of extensive testing of the model, its generality 
can not be established. 

In the present study, only the absorption unit has 
been modelled. However, it is desirable to model the complete 
Benfield system which includes a regenerator and a few heat 
exchangers in addition to the absorber. Since several 
industrial units use promoters other than amines, it may be 
desirable to study their effects on CC> 2 absorption. 
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APPENDIX A 


ESTIMATION OF REACTION RATE constant 


The kinetic constant in dilute , unpromoted solution at 

lew temperature is given by 

log k-j-j = 13.635 - 2895/t (A.l) 

infinite dilution 

which correlates the infinite dilution data very well in the 
range 0° — 40°C, At room temperature, the reaction rate 
constant, is a strong function of ionic strength and is related 
to its value at infinite dilution as 


log kjj = 0.081^ (A. -2 ) 

infinite dilution 

The activation energy is same as that at infinite dilution. 
The follow ing equation for kjj is reasonably well established 


and can be used with confidence at all ionic strengths of 
interest, as well as, in the temperature range of 0 - 110 C, 


log k ZI 


- 13.635 - + 0.08 I ± (A.3) 



